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Heat  Era safer  in  a  Laminary  Boundary  Layer  cm  a  Plata 
with  Consideration  of  Chemical  Nonequilibrium 


by  - 

V.  P.  Stulov 


\ 

fl  *-  *  tT* 


Dearlbed  are  results  of  calculating  the  field  of  flow  and  heat  transfer  in  a 
boundary  layer  on  a  flat  plate  by  the  method*  described  in-l) #  Obtained  were  appro¬ 
ximate  fermlaa,  representing  the  effect  of  variability  of  physical  parameters  on  the 
mafflitude  of  thermal  flov  on  the  wall* 

Equilibrium  Flov*  Heat  transfer  on  a  vail  in  a  boundary  layer  of  reacting  atcnde* 
molecular  mixture  deponds  upon  the  change  in  physical  parameters  ttj>  (product  of  vis¬ 
cosity  and  density  )»  Brand tl  number  P  and  Levis  number  L  across  the  boundary  layer* 
Calculation  of  the  thermal  flov  on  the  vail  in  case  of  equilibrium  dissociation  is 
connected  vith  the  solution  of  a  complex  nonlinear  system  of  ordinary  diffecentlal 
equation s+ 

l 

The  females  gives  below  allow  to  make  an  approximate  calculation  of  the  effect 
of  variability  of  physical  parameters  on  the  thermal  stream  on  the  wall  provided 
these  parameters  change  slightly  in  the  boundary  layer* 

The  expression  for  the  thermal  flow  on  the  wall  with  consideration  of  binary  dif¬ 
fusion  has  the  farm  of 

+ pD  [NO— ■  (1.1)  | 

Changing  over  to  dimensionless  variables*  we  will  present  through  a  dimension**', 
lean  heat  exchange  factor 

.  (1.2)  •' 

VW  "  +  <*w (14) 


Bara  N  -  Buaaelt  number*  R  -  Reynolds  number*  H  -  total  enthaliy*referrad  to  h-e# 
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Solutions  for  h  endcr  ware  obtained  in  farm  of  linear  analyses  by  small  parameters 

I 

ft  characterizing  the  variability  of  physical  parameters  and  the  degree  of  nonequi- 

p  _ 

librium,  lb  king  this  under  consideration,  ve  will  present  the  ratio  in  following 

I 

forms 

(i+SA,)  .  (1.4) 


JL  JL 
rx~v~K 


Haro 


l 


A,= 


_ _  «A. 

“'“1C" 
a  _(«**+ -iKj. 

a*  - “T* - —  '  » 


The  values  &1  represent  the  relative  effect  of  variability  of  physical  parameters 

j 

on  the  heat  flow  on  the  wall*  In  first  approximation  the  nonaquillbrium  in  a  gps  per- 
tiole  does  not  affeot  heat  transfer  (provided,  the  concentration  of  atoms  on  the  wall 
la  equal  to  its  equilibrium  value  at  wall  temperature)* 

The  derivatives  h*iw  are  found  from  solving  the  system  of  first  approximation,  for 
enthalpy*  On  the  basis  of  numerical  calculations  were  obtained  expressions  for  these 
derivative  through  noadlnenalcnal  parameters  of  zero  approximation  (h^,!*,  K  ■  oW 
hoo)  and  values*  characterizing  the  change  in  P,  L  numbers  across  the  boundary  layer* 
is  characteristic  values^x^  and  P  were  selected  magnitudes  of  these  values,  on  the 
wall*  Formulas,  determining  the  effect  of  variability^!!^  and  F  on  the  dimensionless 

heat  exchange  factor  on  the  wall,  will  be  Written.  as  follows  ^ 

»  ,  .  .  ‘‘  *  ■  > 

.  _o.3(^00p0>/|i^.-i). ■  ■■  Oi^/y-i) 

i  *”  “  ^(t-h^+thKVF1)  ..(1^) 

The  affect  of  the  L  number  on  heat  transfer  is  expressed  a a  follows* 

•  A,  =  KP  -£-am(Lw  —  i)+  (*•<) 

The  first  addaat  gives  hare  an  addition  to  the  heat  flew  for  the  ease,  whan  !»■ 

I*  1*  Bda  ease  was  developed  by(if]* 

For  small  Iy-1  formula  frcm[2j  can  be  presented  in  form  of 

A,  =  0.52^— Oe*(A» — 1) 


«fl| 
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Iq  this  way*  for  *  0*6  -  0*7  the  fizet  addend  (1*6)  offers  a  value  by  3°"4W 
higher*  than  formula  (l*7)« 

The  second  addend  (l«6)  is  structurally  analogous  to  (1*5) I  it  takes  into  const- 

I 

deration  the  effect  of  change  in  L  number  fraa  the  value  Iy  across  the  boundary  layer* 

t 

Calculations  were  made  for  two  gases  (ooygen  and  nitrogen),  The  conditions  on  the 

wall  and  in  the  unccndng  flow  changed  in  tho  following  boundaries  t 
I  Sea  attached  3*  ^  , 

I 

i  .  •  .  ; 

Condition  (4)  •  !*••  the  temperature  range*  was  selected  so*  that  in  the  range  of 

•  •  * 

pressures*  according  to  condition  (l)9  the  gas  in  the  stream  was  everyvheres  partially 
dissociated. 

(5)  Conditions  (2)  and  (4)  assure  the  following  boundaries  for  K* 


m 


;  o.i</:<io 

Tor  ranges  of  variables  (l)-(5).  fonoilas  (l*5)  -  (1*6)  represent  satisfactory 

approcdLmations  of  mathematical  calculation  results*  For  a  majority  of  variants  the 

-  .  ♦  • 

discrepancy  docs  not  exceed  10-2 OX,  It  must  be  remembered*  that  the  errors  of  the 
approximation  method  of  solving  the  nonlinear  system  of  aquations  jl^  have  the  vary 
same  order*  In  case  when  temperature  distribution  within  the  boundary  layer  has  a 
maximum*  in  fcrmla  (l^5)  -  (1,6)  instead  of^>  ,P,L  per  infinity  should  be  taken  the 
values  of  these  magnitudes  at  maximum  temperature* 

..  The  first  one  of  the  formulas  (1*5)  is  Is  satisfactory  agreement  with  the  result 
obtained  by|2^J  for  eases*  when  the  fractional  denominator  does  not  exceed  unity*  which 
corresponds  to  small  velocities  of  oncoming  flow* 

The  finding  of  more  accurate  values  of  heat  flow  on  the  wall  with  consideration 
of  variability  of  physical  parameters*  is  limited  by  the  fact*  that  the  transfer  coef¬ 
ficients*  included  in  these  parameters*  are  known  as  functions  of  temperature  and 
concentration  only  vary  approximately*  An  analysis  of  this  problem  is  given  in^* 
According  tofc  the  P  number  for  reacting  atccnic*colacular  mixture  changes  very 
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slightly  with  tbs  change  la  temperature  and  concentration  la  broad  limits* 


■em  -tm 


o  -MOf 


2tax  the  second  f  amnia  is  evident  (1*5)  that  the  correction  to  the  heat  flow  is  in 

i 

,  »  * 

significant*  According  to  calculation  data  *  for  the  selected  dependences  of  F  number 
upon  tenparatur*  this  correction  does  not  exceed  3-5X» 

The  effect  of  L  number  on  the  hast  flow  with  consideration  of  its  variability 
across  the  boundary  layer  reaches  20-30%,  Of  self-sustained  interest  is  a  mors  accurate 
study  the  enthalpy  and.  heat  transfer  profile  in  the  boundary  layer  with  variable  num¬ 
ber  L  at  simplifying  assumptions)^.^  *»  const*-  F  *  const*  equilibrium  dissociation* 
Formla  (1^6)  can  be  considered  as  a  method  of  approximate  evaluation  of  the  imfluenos 
of  tba  L  number  on  heat  transfer* 

Iferi2*  Nonequilibrium  Flow*  The  nonequilibrium  field  was  considered  in  accordance 
with  formulas  frcm^L^*  The  concentration  of  atoms  on  the  wall  was  accepted  as  equal¬ 
ling  its  equilibrium  value  at  wall  temperature*  so  that  the  deviation  from  equilibrium 
appeared  on  account  of  par  tide  deceleration  in  the  boundary  layer  and  diffusion  of 
atoms*  doe  to  concentration  gradients.  It  was  discovered,  that  everywhere s  in  the 
...  boundary  layer  with  exception  of  the  narrow  zona  at  the  wall  (  for  a  sufficiently  low 
degree  of  dissociation  on  the  wall)  there  is  a  lag  in  dissociation*  Examples  of  dig. 

tributing  the  none quillbr ium  share  of  atom  concentration  in  the  boundary  layer  are 
shown  in  fig»l  and  2*  Conditions  on  the  wall  and  in  oncoming  flow  are  listed  in  table- 
la  In  last  column  is  given  tha  scale  by  the  amis  ®f  the  abscissas* 
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^  XAider  the  zone  of  relaxation  for  dissociation  in  boundary  layer  we  will  understand 
such  a  length,  on  which  maxima  deviation  from  equilibria*  dies  out,  e*g.,  up  to  5&* 

I 

Tbs  dependence  of  relaxation  zone  length  lp  upon  outer  conditions  can  be  approximately 

t  • 

presented  in  following  form,  obtained  on  the  basis  of  calculating  a  series  of  variants  t 

1  =  cp-h^,^  (2.1} 

For  the  range  of  temperatures  mentioned  in  por.l^  the  dependence  lp  upon  wall 
temperature  and  sneering  flow  may  be  disregarded* 


Die  linearization  method  to  find  the  none qullibr ium  share  in  the  boundary  layer 
places  a  restriction  on  the  range  of  concentrations  and,  consequently,  an  the  range 
of  temperatures*  The  method  offers  satisfactory  results  at  a  condition  when  0.14=.<£ 
4.0*9  everyvheres  in  the  flow.  The  inequality  id  connected  with  the  fact  that  for 
4^  0.1  and  ct  ^0.9  the  dependence  ^  (T)  has  essentially  a  nonlinear  character* 

The  profile  ctp  acquires  greater  curvature  and  the  last  formula  (1.20)  of  jV]  gives 
unnaturally  greater  values  of  the  nonequilibrium  share  of  atests* 


Assuming  that  a  catalytic  wall  assures  constant  concentration,not  equal  to  its 
equilibrium  value  at  wall  temperature,  i.ei cc^  ^  de(Ty).  The  heat  flow  in  this  ease 
on  the  wall  differs  fraa  its  value  for  equilibrium  dissociation  by  a  magnitude 

A«~  —k ~+w?yW  <*•»= *  (r-  «.  M.  ».  =  * (rw. «.))  (2-2)  r  . 

If  ay 4.  a0(TVr)  (aspacially,  a^O  -  total  recombination)  then  0,  i*e*  heat 
flow  increases* 

m 

is  already  mentioned,  to  calculate  the  effeot  of  none qui librium  in  tha  gas  molecule 
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(partlole)  it  la  neoessary  to  calculate  the  second  approximation  tgr  *  It  ia  appa- 

i 

rest,  that  the  nonequilibriun  in  the  gas  par tiela  has  a  very  insignficant  effect  on 
the  magnitude  of  the  thermal  flow* 

£hr«3*Ncce<iullitoltsa  field  in  the  ftroblen  of  a  Chemically  Insulated  Thermometer* 

The  approximate  dependence  of  the  length  of  the  relaxation  zona  upon  pressure  and  rata 
of  oncoming  flow  in.  the  casa  of  chemically  insulated  thermometer  has  the  fora  of 

I  /p  —  (3-1) 

In  this  ease  lp  rises  with  the  increase  in  u«*>  considerably  slower*  than  in  the  oase 
of  a  catalytic  wall  at  constant  temperature.  Exnaples  of  profiles  of  nonequllibrlin 
chare  of  stems  are  given  in  fig»3.  for  parameter  values  of  oncoming  flew, shown  in  tabic 


i 


111  calculations  were  made  on  the  STRKTA  computer* 


The  author  expresses  thanks  to  C»F«Talenin  for  supervision  of  the  experiment. 

Submitted  July  6,1961 
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